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PACS 67 . 25 . dk - Vortices and turbulence 

PACS 47 . 37 . +q - Hydrodynamic aspects of superfluidity; quantum fluids 
PACS 67. 10. Fj - Quantum statistical theory 

Abstract -Kelvin waves propagating on quantum vortices play a crucial role in the phenomenol- 
ogy of energy dissipation of superfluid turbulence. Previous theoretical studies have consistently 
focused on the zero-temperature limit of the statistical physics of Kelvin-wave turbulence. In this 
letter, we go beyond this athermal limit by introducing a small but finite temperature in the form 
of non-zero mutual friction dissipative force; A situation regularly encountered in actual experi- 
ments of superfluid turbulence. In this case we show that there exists a new typical length-scale 
separating a quasi-inertial range of Kelvin wave turbulence from a far dissipation range. The letter 
culminates with analytical predictions for the energy spectrum of the Kelvin-wave turbulence in 
both of these regimes. 
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Introduction. — Quantum turbulence refers to the 
study of fluids such as 3 He and 4 He that remain liquid 
even at very low temperatures where quantum effects be- 
gin to manifest themselves. In fact, below a characteristic 
temperature denoted as T c and T\ respectively, these flu- 
ids acquire a superfluid component whose viscosity is zero. 
Although T\ is small under usual conditions on Earth, it 
has been discovered recently that superfluids can also be 
found in the core of neutron stars where the temperatures 
reaches values of the order of 10 9 K due to the extremely 
high density. This raises the prospect that this state 
of matter is much more widespread than previously be- 
lieved [1]. Having an inviscid component, one can expect 
that superfluids should easily reach a turbulent state. Un- 
like classical fluids where the size of the turbulent eddies is 
a statistical property that can vary from the small vortices 
in a stirred cup of coffee to hurricanes in stormy weather 
systems, the vortices of turbulent superfluids are tightly 
constrained by the laws of quantum mechanics: the cir- 
culation around a quantum vortex is quantized to integer 
values of the quantum circulation k = 2-kK/M where M 
in 4 He is the mass of the atom and in 3 He is the mass of 
a pair of atoms. It has been established experimentally 
by Vincn [2] that the typical turbulent steady state con- 
sists of a complex tangle of quantum vortices as had been 
speculated by Fcynman [3] and later clarified through the 
numerical simulations of Schwarz [4]. 

As was already understood by Tisza [5] and Landau [6] , 



superfluids can be modeled as consisting of two interpene- 
trating fluids with two different velocities v s (superfluid 
component) and v n (normal component). The relative 
density of each component, (p s and p n respectively), de- 
pends on the temperature of the fluid. In the limit T — > 0, 
p n — > and the fluid becomes totally inviscid. This im- 
mediately raises the question of how a turbulent steady 
state can be achieved at all with energy injection at the 
large scales but no obvious dissipation mechanism at the 
small scales. The current understanding of this issue is 
that the appearance of quantum vortices introduces an 
additional scale (the intervortex mean distance) denoted 
by £. For fluid motions of scales much larger than £ the 
quantization of the vorticity is irrelevant, and the turbu- 
lent spectrum resembles that of a classical fluid with an 
energy flux towards smaller scales. When the energy cas- 
cade reaches scales comparable to I the dynamics of the 
quantized vortex lines becomes relevant. These vortices 
continuously reconnect with each other generating in the 
process Kelvin waves on the vortex lines themselves. For 
scales smaller than £ the interaction between vortex lines 
can be neglected, and then one can study the spectrum 
of Kelvin waves on isolated vortices. Kelvin waves on iso- 
lated vortices are nonlinear, producing an energy cascade 
again towards even smaller scales. When the scale of these 
waves becomes comparable to the core radius of the quan- 
tized vortex, it is believed that elementary excitations (dif- 
ferent in 3 He and 4 He) radiate out resulting in a dissipa- 
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tion mechanism. Due to their crucial role, the properties 
of these Kelvin waves have been studied intensely recently 
in the limit of weak wave turbulence and at T = 0. In par- 
ticular, their energy spectrum at T = was found in [7,8], 
see eq. (7) below. 

So far all of the theoretical studies of the statistical 
properties of the Kelvin waves have entirely focused on 
the limit T = and almost nothing is known for finite 
temperatures. However as soon asT^O there is a friction 
between the normal component and the supcrfluid. This 
so-called mutual friction is parametrised by a scalar a and 
offers another avenue for small scale energy dissipation via 
thermal damping. Having in mind that it is difficult to 
reach temperatures very close to experimentally, gen- 
eralizing the spectrum of eq. (7) for finite temperatures 
represents an essential step towards the formulation of a 
theory of superfluid turbulence. 

Basic equations. — The study of non-linear inter- 
actions of Kelvin waves in the framework of the Biot- 
Savart equation (which determines the velocity field in- 
duced by vortex lines) goes back to the pioneering paper 
of Barenghi, Donnelly and Vinen [9] and has attracted 
a lot of studies since, see [10-14] for different examples. 
Restricting ourselves to wavelengths A <C t, the complex 
and ever-changing topology of the tangle of quantum vor- 
tices can be ignored and we can focus on the behavior 
a single isolated vortex. In this case further progress 
can be achieved by presenting the Biot-Savart equation 
in its Hamiltonian form [15,16]: indw/dt = SH/8w*, 
where w{z,t) = x(z, t) +iy(z,t) denotes deflections away 
from a straight vortex aligned with the z direction, the 
superscript " * " denotes complex conjugation and 
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\J{?\ - Z 2 ) 2 + \w(z!,t) - w{z 2 ,t)\ 



-.dz\dz 2 . (1) 



Here the deflection angle is w'(z,t) = dw/dz. For 
w'(z,t) <C 1 one can then follow a long and cumbersome 
theoretical procedure [7, 12, 17] where the Hamiltonian is 
expanded in power series, represented into the k space, 
w(z,t) — > Wk(t) = Ofe(t)/-y/K, and subjected to a weakly 
non-linear canonical transformation, bk = a k + ..., in order 
to obtain finally an effective (1 3)-wave Hamiltonian 1 

n 4 = \ iK ,2 - 3 blb kl b k2 b k3 + c.c] , (2a) 

(2b) 



fe=fel+fe2+fe3 



1 9 ■> —3* 
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Here the dimensionless parameter \I/ < 1 depends on the 
driving mechanism of the Kelvin waves turbulence and it 
will be considered as an external parameter of the theory. 

The statistical description of Kelvin wave turbulence is 
written in terms of dynamical equations for the Kelvin 



1 Notice that "Hi has the meaning of the interaction energy density 
of Kelvin waves (per unit of vortex length), normalized also by the 
He-fluid density p. Dimensionally [Hi] = cm 4 sec -2 . 



waves action n k = £(\bk\ 2 ) /(2ir) (C is the linear size of 
the periodic box) which are a classical limit of the occupa- 
tion numbers n^: n k — > hn k for rife ^> 1 [18]. Our starting 
point in this paper is the kinetic equation for rife, supple- 
mented by a mutual friction dissipation proportional to a 
temperature dependent dimensionless coefficient a [9] 
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= Stfe {rife/} - aukrik , w fc = A«;fc 2 /(47r) , (3) 



Here u>k is the Kelvin wave frequency of wavenumber k 
and the parameter A is given by A = \n(£/a) correspond- 
ing to the ratio between £ and the vortex core diameter 
a. In typical experiments involving both 3 He and 4 He, A 
varies between 12 and 15. For small temperatures, when 
a <C 1 [19] the mutual friction dissipation auiknk is still 
important being compatible with St k {n k i} <C riktOk for 
weakly interacting waves. Note that the frequency shift 
induced by mutual friction [9] is only a small correction 
with respect to the u>k presented in eq. (3) and can safely 
be ignored at low temperatures. 

Notice that eq. (3) ignores thermal fluctuations of 
Kelvin waves which in thermodynamical equilibrium with 
temperature T are given by the Rayleigh- Jeans distribu- 
tion nbath(^) — T/puik- Here p is the density of liquid 
helium (0.145 g/cm 3 for 4 He) In the framework of stan- 
dard approach [18] and using chain of relations w(z, t) =>■ 
u>k(t) => flfc(i) =>■ bk =>■ nu = "bath(fc) we found that the 
mean squared deviation, caused be thermal excitations, 
(M 2 ) bath i s dominated by waves of length A ~ £, and get 
the estimate of the dimensionless ratio 
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In this estimate for 4 He we took T ~ IK and i ~ 10 -3 cm. 
Since R <C 1 the only relevant mode of excitation of the 
Kelvin waves is via reconnection events between crossing 
vortex lines. 

Finally, for the collision term in eq. (3) with the Hamil- 
tonian (2) was derived in [7]: 

see eq. (5) on the next page 

The collision term in eq. (5) preserves the total energy 
of the system. This allows us to rewrite the kinetic 
eq. (3) as a continuity equation for the energy spectrum 
E{k) = 2u>knk- Integrating the equation over k annuls 
the contribution of the collision term. This means that, 
as usual, there exist a flux of energy down the scales, and 
one can write a continuity equation of the form 2 



rfe(fc) 
dk 



= —au)kE{k) , 



e(fc) = -2 1 Stfe- uj k ,dk' . 
Jk'<k 



(6a) 
(6b) 



2 Considering an isotropic vortex tangle with mean intervortex 
distance £ one relates the energy (density per unit of vortex length, 
normalized also by the density p) flux e ([e] = cm 4 sec — 3 ) over the 
ID vortex line to the bulk 3D fluid energy (density per unite mass) 
flux e ([e] = cm 2 sec — 3 ) with the bridge e = e/l 2 . 
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-3 ( ^- - — - — ) 5 (k 1 -k-k2-k 3 )5(uji - cj k - uj 2 - w 3 ) n fe nin 2 n 3 , nj=n hj , uj^ = u k] . (5) 
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This equation has a stationary scale invariant solution in 
the zero temperature limit when a = 0. This solution 
corresponds to a ^-independent energy flux equal to 
the energy influx eo- It was found in [7] and carefully 
analyzed in [8] : 



E (k) = 



Act, 



1/3 



$2/3 £5/3 



0.304. 



(7) 



Discussions of alternative spectra of Kelvin waves [12, 13] 
can be found in [14] and references therein. 

At finite temperatures eq. (6a) contains terms that de- 
stroys the scale-invariant solution, producing a character- 
istic scale fc*. To find this scale we first use the fact that 
the collision integral converges [8] and make a dimensional 
estimate of de(k)/dk as —e(k)/k. By balancing the LHS 
and the RHS of eq.(6a) using the form of E{k) at T = 0, 
one immediately finds the crossover wave vector k+: 



K = 



3/2 



(8) 



For k <§C fc* the mutual friction term on the RHS of eq. (3) 
is negligible compared to the collision term, and the zero 
temperature solution is expected to hold up to some cor- 
rections that are computed in the next section. We will 
refer to this region of scales as the quasi-inertial range 
of Kelvin wave turbulence. For k 3> k± this solution no 
longer holds, and we need to find a new solution in the 
far-dissipation range as is explained in the last section. 

Energy spectrum in quasi-inertial range of 
Kelvin wave turbulence. — In the quasi-inertial re- 
gion the energy flux e(k) varies slowly with k. This allows 
us to estimate e(k) directly from eq. (7) by substituting 
there e(k) instead of the k- independent eo- This gives the 
following algebraic approximation for e(/c): 



:(fc) = * 2 fc 5 £ 3 (fc)/(C* LN A K ) 3 



(9) 



In the incrtial range of scales this equation is equivalent 
to the solution (7). We note that in general the kinetic 
equation with a = displays another solution which is the 
equilibrium state e(k) = and E k =Const. We can rewrite 
an equation that is dimensionally equivalent to eq. (9) 
that contains that solution as well. This is achieved by 
replacing E 3 (k) by ~^kdE 3 (k)/dk, where the numerical 
factor is chosen to reproduce solution (7) for a = 0. Now 
we have the differential approximation for the energy flux 



e(k) 



dE 3 (k) 



5(C LN A K ) 3 dk 



(10) 



which is a step forward in comparison with the algebraic 
approximation (9), because it accounts not only for the in- 
ertial range solution (7), but also for the thermodynamical 
equilibrium solution E(k) =Const. For k < fc* the mutual 
friction term is small and only slightly modifies the iner- 
tial range solution (7). Therefore we seek a solution of 
eqs. (6a) and (10) in the form: 



E(k) 



Eo(k)f(y), y 



15 / k 

TVfcT 



4/3 



(lla) 



in which /(0) = 1. Note the natural normalization of k by 
the crossover value k+. The exponent and pre- factor in the 
expression for y was chosen such that the final equations 
below [cf. eq. (lib)] for the function f(y) would appear 
as simple as possible. Plugging eqs. (11) into eqs. (6a) 
and (10) we get a second order ordinary differential equa- 
tion for f(y): Vf(y) = 1, where 

Vf(y) ee Ayf(y)f\y) + 8yf'(y) 2 - ll/(y)/'(y) . (lib) 

All the system parameters have dropped out in this equa- 
tion, making it universal. 

We could not find an analytical solution of this equa- 
tion; Instead, we have devised an iterative technique that 
provides approximations to the exact solution with expo- 
nential convergence using a proper choice of a trial func- 
tion. Using the boundary condition /(0) = 1, eq. (lib) 
dictates: /'(0) = —1/11. This suggests a simple form of 
the initial trial function fo{y) = 1 — y/yo with j/o = H- 
The parameter yo can be found in a different way from 
the global requirement 



rvo 

/ f(y)dy = 5 
Jo 



(12) 



This constraint expresses the fact that in the stationary 
case the rate of energy injection is equal to the rate of en- 
ergy dissipation. Simple algebra shows that this condition 
is satisfied with the function fo(y) only if yo = 10. Since 
two different ways to determine yo resulted in yo = 10 or 
11 one is encouraged to improve the accuracy by proper 
polynomial corrections: 



m 

U(y) = (1 ) a m (— ) , a = l 

yu m =0 



(13) 



We iterate this equation starting from the zeroth order 
approximation, fo{y) with yo = 10 which was determined 
from the global requirement (12). Consider next the 2 nd 
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Fig. 1: We present a raw comparison of the absolute value 
of the difference between the LHS and the RHS of eq. (fib) 
for different levels of approximation. The dips correspond to 
points where both sides of the equations match exactly. 



approximation /2(y), in which coefficients a\ and a 2 are 
determined by the requirements 



/'(OH -1/11, f(y ) = -l/. 



(14) 



that directly follows from eq. (lib), while yo again from 
eq. (12). Even though these conditions can be inverted in 
order to exactly determine the coefficients, their precise 
analytical form is not useful here and we prefer to report 
their numerical values in Table 1. One can already see 
the excellent convergence properties of our trial function 
since the coefficients quickly decrease and the value of yo 
is already very close to its naive estimate obtained from 
the zero-order approximation fo{y). 

Finally, we move on to our last level approximation 
fz(y) in which we keep an additional coefficient 0,3 in the 
expansion of the trial function. In this case, one can ex- 
press ai , a2 and as analytic functions of yo that is taken 
now as a free parameter. Its value is obtained by setting 
up an optimization problem in which we minimized the 
mismatch function 



X = 



I/O 



[Vf(y)-l] 2 dy/y 



1/2 



(15) 



The numerical values of the coefficients for all three levels 
of approximation are displayed in Table 1. To illustrate 
the rapid convergence of the trial functions /„ we present 
in Fig. 1 plots of \T)f n (y) — 1|, the absolute error between 
the LHS and the RHS of eq. (lib) for different level of the 
approximations. 

The conclusion of this part is that the energy spectrum 
of Kelvin-waves in the presence of dissipation in the dif- 
ferential approximation (10) can be found with excellent 
accuracy. Moreover, very reasonable accuracy is achieved 
at the first step (with the function /□), that can be written 
in dimensional form as follows: 



E(k) = C L 



A, 



,1/3 



^2/3 £.5/3 



4/3 



(16a) 



Table f : Value of the coefficients for the various levels of ap- 
proximations. 



Approx 


fo 


h 


h 


yo 


10 


9.644 


9.651 







0.123 


0.123 


a 2 





-0.025 


-0.035 


a3 








0.011 


X 


HT 1 


3 x 10~ 3 


5 x 10~ 4 



where the position of the final point ko is given by 



(I) 



3/2 



(16b) 



3/4 



where A = ^327r/3^ ps 14. This predicted spectrum is 
compared to the spectrum at T = in Fig 2. Note that 
this is a powerful prediction since all of the temperature 
dependence via the parameter a is absorbed in ko- This 
means that the energy spectrum, as well as its correspond- 
ing energy flux 



e(k) 



fo 
5 



1 



k \ 4 /3 



/fc\ 4 /3-|2r , k \ 

\~ko~) \ l 5 ~\k~ ) 



(16c) 



can be expressed as universal functions of the ratio k/kg. 

Thanks to eq. (16b) we can estimate ko by relating the 
3D-rate of energy (per unit mass) dissipation e = v'(kL) 2 
to the effective kinematic viscosity v' and to the vortex line 
density L ~ l/£ 2 , see e.g. eq. (1) in [20]. Next, according 
to footnote 2 wc connect s ~ eoj 



kol 



As](v>Ik)^ 



2 and get finally: 
0.03 



a 



3/4 



(17) 



In the "Manchester spin-down" experiment in 4 He [20] it 
was found that for T below 0.5 K (i.e. for T/T\ < 0.25) 
i/~3-10~ 3 k. Taking for concreteness ~ 1 and A ~ 15 
we get the numerical estimate given in eq. (17) and plot- 
ted in the inset of Fig. 2. We conclude that for T ~ IK 
(a ~ 10 -2 ), we have kol — 1 meaning that the Kelvin 
wave cascade disappears and makes way for the possi- 
bility of thermally exited waves. According to the esti- 
mate (17) the cascade reaches the vortex-core scale 1/a 
for very low temperatures T ~ 0.07 K (a ~ 10~ 10 ). For 
even lower temperatures, the energy flux should predomi- 
nantly convert into the production of quasi-particles (such 
as phonons) which may be observed experimentally 

Reaching the edge of the window of locality. — 

It becomes important at this point to remember that the 
original collision integral eq. (5) is convergent under the 
assumption that the wave action spectrum behaves as a 
power law whose exponent is constrained to the following 
interval n(k) oc k~ x ; 2 < x < 9/2. When T = 0, the 
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Fig. 2: The complete wave action spectrum of Kelvin- wave tur- 
bulence at finite temperatures. The symbol £ is a shorthand 
notation for £ = 47rCLjve 1/ ' 3 fc 11/3 /v[' 2/3 . In this representation 
the horizontal line at I corresponds to the Kelvin wave spec- 
trum at T — 0. The gray line in between the crosses is an inter- 
polation between both spectra in the crossover region k « ko- 
The inset corresponds to the estimate (17) for ko- 



exponent x = 11/3 clearly falls within the range of this 
window of locality. However our new corrected spectrum 
in the quasi-inertial range for T ^ is expected to display 
an increase of the local slope (that can be considered as 
a "current" exponent) as one gets closer and closer to fco. 
This suggests that there must exist a fc n0 n-ioc < &o at 
which the wave action spectrum no longer satisfies the 
constraint of locality. This maximal wavevector fc n0 n-ioc 
can be determined by approximating the new wave action 
spectrum as a local power law with a fc-dependent local 
slope which is given by the logarithmic derivative 



lim 

k— >ipko 



d\nn(k) 
91nfc 



11 
~3~ 



1 



4 V 4/3 

IT 1 - ^,4/3 



(18) 



This allows us to demonstrate that locality will be lost 
when x n on-ioc = 9/2 which immediately translates into 
fcnon-ioc ~ fco/2. This result means that our new wave 
action spectrum can only be realized physically up to 
wavevectors about half as the predicted terminal wavevec- 
tor ko- In the next section we show how to go beyond 
this level of description by deriving an extension of the 
wave action spectrum for wavevectors that are much larger 
than fco. The intermediate regime of wavevectors of the 
order of fco is left f° r discussion in the conclusion. 

Energy spectrum in the far-dissipation range. 

The differential approximation was based on the assump- 
tion that the energy spectrum is at least twice differen- 
tiablc function of fc. Nevertheless it lead to the spec- 
trum (16) in which the second derivative does not exists 
at fc = fco- This means that in the vicinity of fco this ap- 
proximation is not self-consistent. This is the reason that 
it leads to the un-physical result that E(k) < for fc > fco. 

In order to find the correct behavior of E(k) in the dis- 
sipative region fc > fco we should return to the full for- 



mulation of the problem (6) with the explicit form (5) of 
the collision integral. To simplify the analysis we analyze 
eq. (5) in the far-dissipative region, where fc 3> fco- In that 
region it is natural to assume that n(k) -C n(fco)- We will 
see that this "bootstrap" assumption is justified since n(k) 
decays exponentially fast. The expansion of St* turns out 
to be a laborious undertaking that requires a degree of 
vigilance; the logical steps are explained in the Appendix. 
We find that the leading term when both these limits are 
taken into account in the form 



-81.6 aC^u^ 



(19) 

where tp is a dummy variable of integration. The collision 
integral must now be balanced by the dissipation due to 
the mutual friction. This leads us to the following equation 
for the wave action spectrum n(fc) 



St* 



-aA«:fc 2 nfe/47r 



(20) 



Restricting ourselves only to decaying solutions and the 
limit of large fc, it is not difficult to see that the general 
solution of eq. (20) can be simplified into 



k 




0.37 fc . 



-^5 ) ' where ( 21a ) 



(21b) 



Note that even though this result is presented here as an 
expansion of the "effective" 4- wave (1 -o- 3) collision inte- 
gral, eq. (5), it does in fact remain identical in the more 
fundamental context of the original 6-wave (3 3) col- 
lision integral of [17]. We are now in possession of 2 dif- 
ferent expressions for n{k). However, since the corrected 
wave action spectrum presented in the previous section 
is only valid for fc <C fco an d the one presented in this 
section is the result of an expansion of the collision inte- 
gral for fc 3> fco, there is no reason to expect that these 
two solutions should be matched in their respective do- 
main of validity. Therefore, we choose the dimensionless 
constant c\ = 1. Nevertheless the crossover region be- 
tween these two different regimes is very narrow around 
the wavevector fc sa fco as one can see from Fig 2. 

Conclusion. — The main message of this letter is that 
the presence of a small but finite temperature is enough 
to dramatically suppress the Kelvin- wave turbulence from 
the simple scale-free eq. (7) to the complex behavior sum- 
marized in Fig. 2. 

Since our theoretical analysis is limited in the crossover 
region, we believe that important additional information 
can be gained from numerical simulations of finite tem- 
perature Kelvin-wave turbulence. However preliminary 
attempts in this direction 3 show that this is a tricky task 



^Preliminary numerical simulations [21] of weak wave turbulence 
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(even for a ID systems like ours) which would require sig- 
nificant efforts in the future. 

The finite temperature results presented in this letter 
appear to constitute the first attempt to bring the theoret- 
ical predictions of Kelvin- wave turbulence closer to realis- 
tic situations. Hopefully these can be used in actual exper- 
iments on supcrfluid turbulence. There are however still 
many questions that lie ahead. For example, how does the 
underlying assumption of weak turbulence, which ignores 
the all-important possibility of quantum vortex reconnec- 
tions, affect the energy spectrum? How to go beyond this 
approximation? 

Appendix: Expansion of collision integral. 

The convergence of the collision integral is ensured only 
through a series of deep cancellations that deserve to be 
explained in more detail. As mentioned in the text, we 
are interested in the limit where one of the wavevector k 
is much larger than the terminal point fcrj. In addition, 
we make the assumption that n(k) <C n(ko). First, let 
us define the small parameter /i = ko/k <C I and for let 
us also refer to the conservation laws of the first line in 
the integral (5) as "spherical" and the conservation laws 
of the second line as "hyperbolic" . In this case, it is clear 
that the hyperbolic geometry requires that k\ be close to k 
while ki and k$ should be much smaller, of the order of kQ. 
On the other hand, the spherical geometry is more permis- 
sive and allows for any of the 3 wavevectors (ki, ki or k$) 
to be close to k. For example, if it is k\ which is of order k 
then it is k% and k% which will be of order fco- In this case 
we can write 

fe 2 = ^fe ; -!<0<I. (A.l) 

Because of the other 2 possibilities, it is now clear that the 
expansion of the part of the collision integral constrained 
by the spherical geometry should be multiplied by a factor 
of 3. Notice that this new prefactor matches exactly the 
prefactor of the hyperbolic geometry albeit with a different 
sign. This feature has already been discussed in [7] and is 
responsible for cancellations of all the terms of order fx 2 , /j, 4 
and /i 6 . Since some of these terms would otherwise be non- 
integrable these cancellations are quite outstanding. For 
the sake of the discussion, let us now go back the spe- 
cial case (spherical geometry) mentioned above in (A.I). 
First of all one should notice that because of the {fc2, — ^2} 
symmetry, all of the odd-ordered terms of the expansion 
identically vanish because of the integration over 0. In- 
terestingly it turns out that all of the terms of order /i 8 
involve quadratic combinations n(k) 2 which must be much 
smaller than n{k) according to our "bootstrap" assump- 
tion stated above. Since the next terms of order /x 9 are 
absent, it is necessary to study terms of order fi 10 . Quite 

reveal that finite size effects play in an important role which is fur- 
ther amplified in this one-dimensional setting. This inhibits wave 
resonances and leads to appearance of several asymptotic ranges. 
This results in great difficulties in choosing the parameters of the 
simulations. 



remarkably, we discovered that there is only one term of 
this order that is linear in n(k); It involves the second 
derivative of the wave action spectrum and is presented 
in eq. (19) 

* * * 
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